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Abstract. The quantum dynamics of a few bosons in a double well potential is studied using a Bose 
Hubbard model. We consider both signs for the on-site interparticle interaction and also investigated the 
situations where they are large and small. Interesting distinctive features are noted for the tunneling 
oscillations of these bosons corresponding to the above scenarios. Further, the sensitivity of the particle 
dynamics to the initial conditions has been studied. It is found that corresponding to an odd number of 
particles, such as three (or five), an initial condition of having unequal number of particles in the wells 
has interesting consequences, which is most discernible when the population difference between the wells 
is unity. 

PACS. 0 3.75.Lm, 05.45.Yv, 03.75.-b 


1 Introduction 

Ever since the ultracold bosonic atoms could be experi¬ 
mentally prepared by confining in a single quantum state, 
called a BEC [HI], the researchers in the field have in¬ 
tensely tried to broaden their research on cold atoms. Ma¬ 
nipulation of the atomic gas to emulate various quantum 
many body phenomena has emerged into an exciting re¬ 
search endeavor of the atomic physics and the condensed 
matter physics community which is aptly complimented 
by the discovery of optical lattices with precisely tunable 
interaction potentials via manipulating the laser parame¬ 
ters and Feshbach resonance. Such effects are enormously 
facilitated by an ultra-clean, phonon-free system, and fi¬ 
nally, a vastly magnified version of the crystal lattice. 

Tunneling of particles through a classically impenetra¬ 
ble barrier is a classic problem of quantum physics [3]. 
When interaction between these particles are included, it 
may help or hinder the tunneling phenomena. The time 
resolved tunneling probability may demonstrate interest¬ 
ing effects of the roles of the interaction parameter(s) and 
the initial configuration of the particles. For example, the 
phenomena of time evolved pair tunneling of the particles, 
as opposed to the individual tunneling across the barrier, 
can crucially depend upon the initial state of the system. 

In the regime where the atoms are weakly interacting, 
tunneling phenomena of individual particles dominate, as 
it is the case for normal Josephson junctions. However 
as the (repulsive) interaction grows stronger between the 
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atoms, two of them located at one side of the barrier can¬ 
not tunnel independently and thus a pair tunneling be¬ 
comes inevitable. There can also be a ‘conditional tunnel¬ 
ing regime’ where tunneling of a single particle can happen 
only in the presence of a second particle that acts as switch 

a- 

A simplified version to study correlated particle dy¬ 
namics in presence of confining potential is to consider a 
single particle or a few particles in a double well potential. 
Albeit straightforward, it has the potential to demonstrate 
a range of fundamental quantum phenomena with regard 
to the tunneling dynamics of the particles and abilities to 
manipulate them in terms of suppression of the tunnel¬ 
ing probabilities, and thereby trapping them in one of the 
wells. Such trapping phenomena are experimentally real¬ 
ized with a BEC The studies involving a few bosons 
have turned out to be more relevant in recent times af¬ 
ter experimental successes of the ‘Boson Sampling’ [7118119] 
where a small number of bosons were used for experimen¬ 
tal demonstration of achieving unprecedented control of 
multiphoton interferences in large interferometers. These 
techniques offer huge prospects of simplifying the quan¬ 
tum computation problem and speeding it up further. 

A two mode approximation, valid when the energy dif¬ 
ference between the two lowest single particle eigenstates 
is far smaller than all other energy states, can describe the 
tunneling between different Bloch bands in an optical lat¬ 
tice nnj. In this work, we shall investigate a two-site Bose 
Hubbard model (BHM), which is the simplest candidate 
to investigate the dynamics of correlated bosonic atoms 
in a double well potential or a bosonic 
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junction m- At the outset, it is helpful to mention that 
we shall mainly focus on the physics of weak and strong 
inter-particle repulsion limits and the sensitivity of initial 
conditions on the tunneling dynamics. 

The dynamical evolution of the Fock space for a simple 
two-site Bose Hubbard model (BHM) (without the den¬ 
sity exchange term) with two bosons has been investigated 
and the tunneling probabilities are computed as a func¬ 
tion of time [T 7 III 8 ) . However a detailed analysis of the 
quantum dynamics in strong and weak coupling regimes 
and the sensitivity of the time evolved state to a variety 
of initial states were lacking. This is particularly relevant 
for engineered waveguide lattices to achieve a certain pre¬ 
ferred final state. Thus it is interesting to consider a few 
(N) bosons {N > 2, for example, = 3,4 etc) as the com¬ 
plexity of dynamics (compared to iV = 2) is inevitable as 
the many body effects will become more conspicuous for 
a larger assembly of particles. Furthermore, an inter-site 
density exchange term in BHM, relevant for a gas of dipo¬ 
lar bosonic atoms, can also be considered in the present 
context. Inclusion of this term in the BHM is known to 
have density ordering effects and is responsible for a rich 
phase diagram in a lattice [T9]. 

Motivated by the above prospects, we have considered 
a few bosons in a Bose Hubbard model (BHM) and inves¬ 
tigated the quantum dynamics. In particular, we have in¬ 
vestigated the dependencies of the tunneling probabilities 
in the strong and weak on-site interaction limits, briefly 
the effect of the density ordering term therein and the sen¬ 
sitivity of the dynamics to a variety of initial conditions. 
Among other results, the time evolved dynamics is seen to 
be crucially dependent on the initial state of the system in 
which it is prepared, particularly when the initial popula¬ 
tion difference between the two wells is unity for an odd 
number of particles. Further, we inlcude a brief discussion 
on the effect of using an admixture of initial states on the 
tunneling oscillations. 

In the following, the presentation of the paper is or¬ 
ganized as follows. The next section deals with studying 
the quantum dynamics exactly for a system consisting of 
a few bosons conhned in a double well potential and de¬ 
scribed by a BHM on a two site lattice. Hence, we present 
our results on the effect of different initial conditions on 
the tunneling dynamics. In particular, we have included a 
brief discussion on using different admixture of states as 
initial conditions. The implications of our results on the 
experimental scenario is presented thereafter. 


2 The Bose-Hubbard Model and the 
tunneling dynamics for a few bosons 

Even though we are going to restrict ourselves to the usual 
(short ranged) Bose Hubbard model [TTII^ , we include an 
extended density ordering term while deriving the equa¬ 
tions of motion (EOM) with the motivation of investigat¬ 
ing its competing effects with the on-site term on the tun¬ 
neling dynamics. It is relevant to mention that for dipo¬ 
lar bosons, such extended range interaction potentials are 


important to include, as the research of ultracold dipolar 
gases gained interest with the experimental realization of 
Bose condensed Cr atoms which hosts large long range in¬ 
teractions m- As will be immediately clear, the extended 
term, unlike that for a lattice, only renormalizes the on¬ 
site interaction for a double well. 

For a system of N interacting bosons occupying the 
weakly coupled low lying energy states of a symmetric 
double well potential, the BHM Hamiltonian is written 
as, 

Hbhm = -J{a\a2 + a^ai) -I- + a^^al) 

-[-V(a\aia\a2) = Hj -\- Hjj ~\- iTy, (1) 

where a|(a 2 ) are the creation (annihilation) operators 
of bosons in the left (right) wells, J being the tunneling pa¬ 
rameter (J > 0) between the two modes, U is the strength 
of the on-site interaction {U > 0) and V is the strength of 
the extended density interaction (or exchange interaction) 
that has, as mentioned earlier, implications in formation 
of density order phases and are suitable in the context of 
dipolar bosons. It may be noted here that all the energy 
scales including the time evolution are expressed in units 
of tunneling frequency, J. It may be noted that the Hy 
term includes Hu in the following way m, 


Hv = V 


-N- 


u 


( 2 ) 


where N = ni +n 2 = ajoi -I- 0 ^ 02 . In the light of this, the 
density exchange term (Hy) has no role, apart from nor¬ 
malizing the interaction strength U to U' (= U — V). This 
consequently implies a noninteracting scenario for U = V. 
The role of the inter-particle interaction becomes only rel¬ 
evant for U ^ V. Further, the negative sign in the expres¬ 
sion of U' means that it can either be negative (attractive) 
for U < V [23] or positive (repulsive) for U > V, for which 
we have considered V = 2.U and V = 0.5C/ respectively 
as the representative values. However due to the symmet¬ 
ric nature of the model, a sign change in interaction does 
not contribute in its dynamics. Albeit, in our analysis we 
have used different values corresponding to attractive and 
repulsive region (as a pathological case) which effectively 
provides two different interaction strengths manifesting 
solely the role of interaction magnitude instead of its char¬ 
acter. Further, we have distinguished the weak and strong 
coupling regimes by assuming U — 0.1 and U = 12 (both 
in units of the tunneling frequency, J) respectively. The 
‘strong’ and ‘weak’ will carry these values along through¬ 
out the manuscript. Other representative values have been 
assumed for the computation of tunneling dynamics, how¬ 
ever they yield no new qualitative inference. 

To obtain the tunneling dynamics, the state vector of 
the system is expanded in the basis of Fock states for a 
constant particle number N, as in the following. 


N 


\m) = E 


Ci{t) 






|o), 


(3) 
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where out of N particles, I are in the left well and N — l 
are in the right well and q {t)s are the complex coefficients. 
The EOM can be written as, 

^ j. (4) 

at 

In terms of the coefficients c/(t), EOM is expressed as, 


i = -kci+i - ki-iQ-i + aici + bid, (5) 

where h = J^{1 + 1){N - 1), ai = + {N - If - N], 

bi = V[l{N-l)] 

For the case two bosons {N = 2), Eq.Q reduces to 
three coupled equations as in the following, 


i 


d 

dt 



/ U -V2J 0 \ /co\ 

-V2J V -V2J Cl . (6) 

V 0 -V2J U / \C2/ 


one in the left, can be denoted by, 
. dci 


= —VNJco — \/2(7V — 1) Jc2 


N'^-m+ 2 


Uci + {N-l)Vci. 


For an equal distribution of the bosons with 7V/2 in each 
well {N: even) is 

= -^VNiN + 2 ) {cN/2+l + CN/2-l) 

+ —N{N — 2)civ/2 + -^N‘^cn/2- 

However for discussing our results in the following section, 
we restrict ourselves to the case of a few bosons. 


3 Physical Observables and Results 


Hence the particle occupation probabilities can be ob¬ 
tained by solving these coupled equations. The various 
initial conditions that can be thought of are, 

co(0) = l,ci(0) =C2(0) =0; ci(0) = 1, co(0) = C2(0) = 0; 

C2(0) = l,co(0) =ci(0) =0. (7) 

In short we shall denote them as (100), (010) and (001) re¬ 
spectively, where (100) means that, initially all the bosons 
are in the right well with the left one being empty. Sim¬ 
ilarly, (010) denotes one in each well, while (001) implies 
both in the left well with the right well being empty. 

Similarly, for the case of three bosons, one gets four 
coupled equations, which are. 


As emphasized earlier, we are interested in studying the 
quantum tunneling dynamics of a few bosons in a double 
well potential. We are mainly interested in the effect of 
the inter-particle repulsion and density exchange for two, 
three, four and five bosons and the variation in the tunnel¬ 
ing dynamics associated with different initial conditions. 
In the following, we describe the cases corresponding to 
two, three and four bosons separately. In this regard, an 
useful (and experimentally measurable) quantity to study 
the tunneling dynamics of bosons can be the population 
in one of the wells (say the right well), PR{t) as discussed 
in the following discussion pO] . 

For a system of two bosons, the right well population, 
PR{t) can be expressed as [IT], 


/co\ 

/ 

31/ 

-V3J 

0 

^ \ 

cr _ 


-y/3J 

to 

-2J 

0 

C2 - 


0 

-2J 

U + 2V 

-V3J 

Vca/ 

1 

0 

0 

-Vsj 

317 / 



( 8 ) 


As earlier, these equations can be solved for four differ¬ 
ent initial conditions, namely, (1000), (0100), (0010) and 
(0001), with implications as before. For example, (0100) 
denotes a situation where two particles are in the right 
well, with the other in the left well and so on. 

A straightforward extension yields similar set of equa¬ 
tions (now a set of five) and the corresponding initial con¬ 
ditions for iV = 4. We have repeated the procedure till 
= 16. For brevity, we skip them here. 

An extension of our results to the case of N bosons is 
possible via the method of induction. A straightforward 
application of this method on Eqs. and m yileds, 

‘A eo, (9) 

where all the bosons are in the right well. Similarly the 
case corresponding to {N— 1) bosons in the right well and 


PH(t) = |co(t)P + i|ci(f)p, (10) 

which is a superposition of the probabilities of both the 
bosons in the right well and half of that corresponding to 
one in each well. 

It may be noted that for U' = 0 {i.e. U = V), Rabi 
oscillation of the particles between the wells is observed 
(Fig. 1(a)). At small values of U', the atoms can still tun¬ 
nel independently, similar to that of the normal Josephson 
junction, however the time period for oscillation becomes 
enormously large which signals the onset of a trapping 
scenario. A similar scenario has been reported by Zollner 
et al. EO], where they have found that for g = 1.3 {g 
being the strength of the pairwise potential), the time pe¬ 
riod is as large as 2 x 10^ s. In this situation, two or 
more atoms residing in one of the wells, form a ‘repul¬ 
sively bound pair’ [24] and hence tunnel together. Such 
phenomena are difficult to contemplate in crystal lattices 
owing to relatively much shorter life times associated with 
the decay processes. 

For the weakly interacting case U = 0.1, it is impor¬ 
tant to note that PR{t) collapses, as evident from Fig. 1(b). 
However there is again a ‘revival’ as time progresses and 
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Fig. 1. (Color online) -Pfi(t), defined in text, as a function of time for two bosons is shown for non interacting case {V = U) 
in (a). In (b) left triangles with dots (black) denote -Pfl(t) corresponding to U' = 0.05 and solid line (red) denotes Pnit) 
corresponding to U' = —0.1. Similarly in (c) dashed lines (red) denote Pnit) corresponding to U' = Q and solid line (black) 
denotes Pnit) corresponding to U' = —12 The time t in the a:-axis is measured in units of the tunneling frequency, J and is 
true for all plots. 


this phenomenon is repeated with increasing time. At time, 
t = 0 , the system is prepared in a definite state (as de¬ 
scribed by the initial conditions in the preceding section) 
and the two terms in Eq. (US are correlated. However as 
time increases, the oscillations corresponding to different 
initial excitations pick up different frequencies and hence 
become uncorrelated, thereby leading to a collapse. With 
further increase in time, the correlation is partly (depend¬ 
ing on the value of U) restored and revival occurs. This 
behavior repeats itself and thus an infinite sequence of 
collapse and revivals are obtained [5^. At large values of 
U, namely, U = 12, termed as the fermionization limit 
(where bosons avoid each other and thus obey an ‘exclu¬ 
sion principle’), there are faster oscillations with smaller 
amplitudes, however PR[t) becomes zero eventually (see 
Fig. 1(c)), signaling a tunneling of the atoms at large time 
scales. The time period of such ‘eventual tunneling’ phe¬ 
nomena, Tr (say) increases with the increase in the inter¬ 
particle repulsion, U, thereby signaling intense trapping 
effects. The analytic expressions for Tr corresponding to 
the noninteracting ({7 = V), and interacting (considering 
two pathological cases V = 2.U and V = 0.5U) cases are 
obtained as. 


Tr = n/J for H = [7 
Stt 

U - \/64 J2 + c/2 
47r 

u - Vie J2 -b c/2 


( 11 ) 

for V = 0.5U 
for V = 2U 


Thus, as a function of U, Tr scales linearly and the agree¬ 
ment between the analytic expressions and the correspond¬ 
ing numeric estimates are shown in Fig. (2). So Tr is in¬ 
sensitive to U values for the non-interacting case, while 
it sharply increases as U is increased for the interacting 
cases. 

Let us now analyze the impact of different initial con¬ 
ditions on the tunneling dynamics of two particles in a 
double well potential. To prepare an initial state with a 
population imbalance, a tilt in the form of a linear po¬ 
tential —r]X {f] > 0) can be superimposed For a 


reasonably large rj (magnitude of the tilt), all the particles 
can be made to reside in one well. The subsequent dynam¬ 
ics can be studied by allowing 77 —>■ 0 within some charac¬ 
teristic time scale. Motivated by such prospects of experi¬ 
mentally creating different initial states m, we study the 
tunneling dynamics subject to different initial conditions. 

In the weak coupling regime {U = 0.1) and the so 
called attractive limit {U' < 0 or F = 2.U), PR{t) oscil¬ 
lates and slowly dampens for the initial conditions given 
by (100) and (001) (see Fig.3(a)), while the damping is 
faster with further weakening of the interaction field (U' > 
0 or F = 0.517) (Fig.3(b)). The situation in the strong cou¬ 
pling limit (U = 12) show Rabi oscillations with different 
frequencies corresponding to the F = 2U (Fig.3(c)) and 
F = 0.517 (Fig.3(d)) situations. It may also be noted that 
if we start with an initial condition ( 010 ) where one boson 
resides in each well, Eq.(3) becomes, 

|!f(l)) = ci(t)aja|| 0 ) ( 12 ) 

It can be shown that the above state is an eigenstate of 
the Hamiltonian (Eq.(l)) and hence the dynamics is frozen 
which is seen from Fig.(3), where PR{t) stays at 0.5 irre¬ 
spective of the values of parameters used. The frequency 
and time period of these oscillations depend upon the in¬ 
teraction parameters used in this work. 

Let us now concentrate on the case of three bosons 
{N = 3). The PR{t) in this case is defined as, 

P^(t) = \co{t)f + licit) f + l\C2it)f. (13) 

Similar to the case of two bosons, here PRit) is a com¬ 
bination of probabilities of all of them in the right well 
(|co(l)P), two in the right and one in the left with an am¬ 
plitude I and one in the right and two in the left with 
an amplitude ^ respectively. While qualitatively the tun¬ 
neling behavior remains unaltered as compared to two 
bosons, with regard to Rabi oscillations at C/' = 0 (not 
shown here) and there is a temporary decay of the ampli¬ 
tude of oscillations (Fig.4(a)) due to the beating phenom¬ 
ena for 17' 7 /: 0. It can also be seen that F = 0.5C/ registers 
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Fig. 2. (Color online) Time period Tr of two bosons as a function of U is shown. In the figure, the circles with dotted line 
(blue), triangles with dotted line (black) and squares with dotted line (red) denote the analytical results for the non-interacting 
case {V ~U), and interacting cases with V = 0.517 and V = 2U respectively. The associated open squares are the numerically 
obtained results. 





(a) (fe) (c) (d) 

Fig. 3. (Color online) Pnit) as a function of time corresponding to three different initial conditions is shown. In (a) U' — 0.05, 
and (b) U' = —0.1 solid line (black) denotes (100), dashed line (red) denotes (010) and diamonds with dots (blue) denote 
to(OOl). In (c) U' = 6 and (d) U' = —12 solid line (black) denotes (100), dashed-dotted line (red) denotes (010), and dashed 
lines (blue) denote (001). See text for details on notations. 


more significant decay of the amplitude owing to trapping 
effects. At large U, (Fig.4(6)), the time period of ‘eventual’ 
oscillations becomes very large. The time period is about 
an order of magnitude larger compared to two bosons. In 
Fig.4(b), for U' = 6, one can observe tunneling phenom¬ 
ena at larger time scales, while at U' = —12 (large V), the 
tunneling of atoms take a very long time and we do not 
observe any tunneling until t = 60 and even to much large 
values of time (not shown here). This indicates emergence 
of trapping phenomena for large values of the exchange 
intercation V in the large U regime. 

A close scrutiny of different initial conditions for three 
bosons in the small U regime reveals an interesting ob¬ 
servation. Pfi(t), corresponding to the initial condition 
(0100) ((0010)), starts with | (i), as expected. However, 
as time progresses, Pfi(t) modulates between values ap¬ 
proximately 0.85 and 0.15 (Fig.4(c) and (d)). Thus in the 
weak coupling regime {U ~ 0.1), the fraction of the to¬ 
tal number of bosons occupying the right well is becom¬ 
ing larger (smaller) than | (|), thereby indicating a ten¬ 
dency of accumulation of particles in one of the wells. This 
seems like an interesting result as an accumulation of par¬ 
ticles is not expected for t/ < J {U = 0.1 in units of J 
here). Pnit) oscillating between values such as, | and ^ 


may have been more commonly expected. The other initial 
condition, namely, (1000) (or (0001)) does not exhibit any 
noteworthy feature and hence not included for discussion. 

Therefore, it indicates that in case of odd number of 
particles, when both the wells contain unequal number of 
particles and the population difference between the wells 
differs by unity, then such a scenario of accumulation of 
particles may be observed. However in the noninteract¬ 
ing limit (with V = U), such accumulation of particles 
vanishes and Pnit) oscillates between | and Similar 
result emerges for large U limit {U = 12) where the ac¬ 
cumulation of particles ceases to be a possibility owing to 
trapping effects. We have skipped these plots for brevity. 

As an extension of the ongoing discussion, we take a 
look at the case of four bosons. Here PR{t) is defined as, 

Pnit) = |co(t)P + J|ci(f)p + i|c2(t)P + i|c3(t)P, (14) 

There is no qualitative difference in the behavior for PR^t) 
both in H = 2U and V = 0.51/ cases corresponding to the 
weak coupling regime between this and those for two or 
three bosons. In the strong coupling case as expected, the 
localization is strong, and a complete tunneling of all the 
particles is prohibited over a very large time scales. Thus 

















































6 


S. Dutta et. al.: Tunneling dynamics of correlated bosons in a double well potential 





(a) (6) (c) (d) 

Fig. 4. (Color online) Pnit) with the initial condition (1000) for three bosons. In (a) square with dots (black) denote U' = 0.05 
and solid line (red) denotes U' = —0.1 and (b) solid line (black) denotes U' — Q and dotted line (red) denotes U' = —12. Pnit) 
with (0100) initial condition is depicted in (c) with U' — 0.05 and (d) with U' = —0.1. 


the notion of (Rabi) oscillations at large U as inferred 
earlier, is no longer observed, at least for time scales ~ 10 ^ 
(in limits of the tunneling frequency). 

Again an inspection of Pnit) with different initial con¬ 
ditions such as ( 10000 ), ( 00001 ) and ( 00100 ) yield results 
similar to those corresponding to ( 100 ), ( 001 ) and ( 010 ), 
respectively, for two bosons with the last one in either case 
yields, Pnit) = 0.5 for all t and this value is fairly insen¬ 
sitive to the values of U and V. For (01000) and (00010) 
as possible initial conditions, that is, three particle in the 
right well and one in the left, Ffi(t) oscillate between val¬ 
ues I and and vice versa as expected. 

Thus the effect of initial conditions seem to be impor¬ 
tant for an odd number of bosons in a double well, spe¬ 
cially when the population difference between the well is 
unity. The claim is substantiated by looking at the case of 
five bosons, for which, at small values of U, three particles 
in the right well and two in the left (or vice versa) produces 
probabilities nearly | and ^ as time progresses, which are 
greater than | and |, thereby indicating a possibility of 
accumulation of particles. However, four particles in one 
well and one in the other demonstrates no such accumu¬ 
lation tendencies, where PR{t) values oscillate between | 
and I as expected. The plots are skipped here for brevity. 

4 Admixture of states 

Further emphasis on the effect of initial conditions can be 
given as follows. Instead of choosing a particular initial 
state, one can consider an admixture of states. For exam¬ 
ple, for the case of two particles, instead of assigning an 
initial state ( 100 ), we may consider an admixture of the 
form, 

^(100)-f 7(010)-b (5(001) 

with the restriction, |/3p -b jyp -b |(5p = 1. In this spirit, 
we have considered small deviations from the pure states 
by suitably choosing /3, 7 and 6 and looked at the time 
evolved states via Pnit) for a comparison with those for 
the pure states. The corresponding plots for some specific 
choices of /3, 7 and 5 are presented in Fig.(5). A slight 
deviation from a pure state, say ( 100 ) as an initial state, 
results in a slightly different (Figs. 5 (a) and (c)) or com¬ 
pletely different dynamics (Figs.5 (b) and (d)) where in 


the latter case, a small mixing of the probability ampli¬ 
tudes corresponding to ( 010 ) state yields an oscillatory 
dynamics. Thus the probability amplitudes of the initial 
states can slightly be modified to yield a desired oscilla¬ 
tory dynamics. 

5 Time averaged dynamics - extrapolation to 
large N 

In order to draw relevance of these results to the experi¬ 
ments done on cold atoms, we need to extend the studies 
for a large number of bosons. An exact computation of the 
quantum dynamics for such a large system is difficult. So 
we present a time averaged Pnit), denoted by a, which is 
defined as, 

a = ^ ^ Pnit) (15) 

by computing the EOM exactly for upto 16 bosons and 
plotted a as a function of 1/A^ {N being the number of 
bosons) and hence the results are extrapolated to A^ —>■ 00 
(or 1/N —>■ 0). Here T is taken as 30 in units of 1/J (J be¬ 
ing the tunneling amplitude). The results in the weak and 
strong coupling limits are presented in Fig.( 6 ). In Fig. 6 (a) 
which corresponds to a small U regime, a does not de¬ 
pend on N and stays at an average value of 0.5, regardless 
of whether interaction effects have been included. In the 
large U regime (in fig. Fig. 6 (b) ), a , although flat, yet 
different for different values of V corresponding to smaller 
number of particles, shows a linear fall off as N becomes 
large. In the limit N —^ 00 , a becomes small 0.1 — 0.2 
(the extrapolated value), re emphasizing the onset of the 
trapping effects as the time averaged probability for the 
particles to spend in one of the wells (right well here) be¬ 
comes low. Hence there is indeed a depreciation in the 
value of the time averaged right well population in pres¬ 
ence of a large number of bosons in a double well potential, 
however qualitatively similar physics can be expected as 
that for a few bosons. Summarizing the above discussion, 
we conclude by saying that the tunneling period not only 
increases with the particle number, but also depends on 
the interaction strength. Saturation behaviour of the time 
period (that is intense trapping) is expected in the limit of 
N for the strongly interacting regime (which is clear from 
Fig.( 6 )). 
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Fig. 5. (Color online) -Pfl(t) with different admixtures of initial conditions are shown for two bosons. In (a) U' = —0.1, squares 
with dashed lines (black) denote the initial condition (100) and solid line (red) denotes the initial condition (VO.9, 0, VO.l), (6) 
U' = —0.1, dashed lines (black) denote the initial condition (010) and solid line (red) denotes the initial condition (VO.l, \/0.9,0), 
(c) U' = 6, squares with dashed lines (black) denote the initial condition (100) and solid line (red) denotes the initial condition 
(v/0.9, 0, VO.l), (d) [/' = 6 here dashed line (black) denotes the initial condition (010) and solid line (red) denotes the initial 
condition (%/0.1, %/0.9,0) A small admixture leads to an oscillatory dynamics in (b) and (d). 



(a) (6) 


Fig. 6. (Color online) Time integrated value of F’fl(t) i.e. q as a function of inverse of the number of bosons is shown for (a) 
weak and (b) strong coupling cases. In (a) circles with dotted line (green) denotes a for V = U (non-interacting),solid line (red) 
denotes a for V = 0.517 and diamonds with dashed line denotes a for V = 2U. In (6) circles with dotted line (green) denotes a 
for V = U (non-interacting), right triangles with solid line (red) denotes a for V = 0.517 and diamonds with dashed line (blue) 
denotes a for V — 2U. 


6 Conclusions 

We have carried out a detailed enumeration, though by 
no means exhaustive, of the effects of onsite inter-particle 
repulsion on the tunneling dynamics of a few bosons in a 
double well potential. The strong and weak coupling limits 
are compared and contrasted with regard to the study of 
tunneling dynamics. Further, the sensitivity of the particle 
dynamics to different initial conditions is closely scruti¬ 
nized. For an odd number of particles in the limit of weak 
repulsion, a population difference of one particle among 
the two wells seems to demonstrate accumulation tenden¬ 
cies. However, no such behavior is observed for the popu¬ 
lation difference to be larger than one. Also the effect of 
an admixture of initial states on the tunneling oscillations 
has been studied. It is premature to comment on the im¬ 
plication of this result to more elegant phenomena, such 


as using it as an ^atom switch' etc, however we feel that 
our results can motivate further experiments in the study 
of atomic dynamics in presence of correlation effects. 
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